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Abstract 

We study the critical behavior of the D {2 < D < A) dimensional Gross-Neveu 
model with a Thirring interaction, where a vector-vector type four-fermi interaction 
is on equal terms with a scalar-scalar type one. 

By using inversion method up to the next-to-leading order of 1/N expansion, we 
construct a gauge invariant effective potential. We show the existence of the chiral 
order phase transition, and determine explicitly the critical surface. It is observed 
that the critical behavior is mainly controlled by the Gross-Neveu coupling g. The 
critical surface can be divided into two parts by the surface 5 = 1 which is the 
critical coupling in the Gross-Neveu model at the 1/iV next-to-leading order, and 
the form of the critical surface is drastically change at 5 = 1. Comparison with the 
Schwinger-Dyson(SD) equation is also discussed. Our result is almost the same as 
that derived in the SD equation. Especially, in the case of pure Gross-Neveu model, 
we succeed in deriving exactly the same critical line as the one derived in the SD 
equation. 


1 Introduction 

One of the most important subject in the elementary particle physics is how to treat 
non-perturbative effect and how to obtain non-perturbative information. Since the work of 
Nambu-Jona-Lasinio(NJL)|l|, chiral symmetry breaking have been studied extensively in 
various models, e.g. QCD, strong coupling QED|^, gauged NJL[^ Q, etc.. Particularly, 
the mechanism of dynamical mass generation has been applied to models beyond the 
^e-mail: dateki@eken.phys.nagoya-u.ac.jp 
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standard model, for example, walking technicolorj^, top mode standard model[^, etc.. 
In these models, scalar 4-fermion interaction plays an important role. Snch a 4-fermion 
interaction in D{2 < D < 4) dimensions, the Gross-Neven model0], has been extensively 
stndied. 

Recently dynamical symmetry breaking and dynamical mass generation in the Thirring 
model[^, which has vector-vector type 4-fermion interaction in D dimensions(2 < D < 4), 
has been stndied by many anthors 0.53.0- Usnally the Thirring model is stndied in 
the form where the fermion has a minimal conpling with a massive vector boson throngh 
introdnction of the vector anxiliary held. Althongh this theory (as well as original Thirring 
model) has no local gange symmetry, several anthors treated this massive gange held as a 
gennine gange heldlQ. Itoh et ah have reformnlated Thirring model as a gange theory 
by nsing the hidden local symmetry ||12|| , and emphasized the advantage of maintaining 
snch a manifest gange symmetry in the analysis of Thirring model in terms of SD eqnation. 
Kondof^ fnrther analyzed snch a formnlation in terms of the inversion method [|I^. 

It wonld be interesting to consider the critical behavior of the model in which the 


above two types of 4-fermi interaction coexist. Actnally, Kim et al.[^ stndied snch a 
model bnt withont manifest gange symmetry. Due to each of the gauge symmetry their 
analysis of the SD equation unfortunately depends on their assumptions. 

In this paper, we analyze Gross-Neven model with Thirring interaction by introducing 
hidden local symmetry and using the inversion methodinstead of SD equation. We 
obtained the generating functional of the next-to-leading order and the effective potential 
by the inversion method. The critical line is obtained in the parameter space {g,gy,N) 
where g and g.^ are the dimensionless couplings of Gross-Neven and Thirring model, 
respectively, and N is the number of fermions. The phase structure drastically changes 
at g=l which is the critical coupling of the pure Gross-Neven model at the leading order 
of 1/A^ expansion. For g > 1 the phase structure is essentially determined by the Gross- 
Neven interaction being only operative in the 1/A^-subleading corrections, while for g < 1, 
the the symmetry breaking take place only due to the subleading effects where the Thirring 
interaction yields attractive force to form the condensate in contrast to the Gross-Neven 
interaction which is opposite. 
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Following the procedure of Ref.lM, we reformulate the original theory as a gauge 
theory in section |^. In section ^ and section ^ we calculate the effective potential for 
the chiral order parameter {'ip'ip), by using the inversion method in 1/A^ expansion. Based 
on this potential, we study spontaneous chiral symmetry breaking and explicitly obtain 
the critical surface which is the boundary between symmetric phase and spontaneously 
broken phase. Comparison with the SD equation is discussed in section |^. 


2 Hidden Local Symmetry 

Let us consider D-dimensional Gross-Neveu model with Thirring interaction {2 < D < 
4). The Lagrangian is given by 

C = + ^(T(a;)T(a:))2 - (1) 

where is a 4-component Dirac spinor with index a {a = 1, 2, ■ • •, N) and belongs to 
a fundamental representation of U{N), and 7 ^(/i = 0,1, 2, • • •, D — 1) are 4x4 matrices 
satisfying the Clifford algebra = 2(y'^,yl. Index a will be omitted in the following. 

Here we introduced the source J = constant in order to study (TT). Using the auxiliary 
held, we may rewrite this Lagrangian in to the form 


= +J)T- 


V- 


1 T T 1 1 

+ -—Al + 


2 G 5 y/N 2Gy ^ y/N 




( 2 ) 


where the auxiliary massive vector held corresponding to the current is not 

really a gauge boson as it stands, and this Lagrangian actually has no gauge symmetry. 

Itoh et aL|T^ introduced hidden local symmetry to the Thirring model and treated it 
really as a gauge theory. Following the same procedure as in Ref. [|^ , we introduce hidden 
local symmetry to our model. What we should do is only the following substitution: 

[ 4/ = 


The Lagrangian d^) is written as 

1 o 


C" = i/j{i ^ + J)i/j 


2G. 


-a 


A^ = A^- VNd^<P . 




( 3 ) 


Vn 


2G 


V 


Vn' 


( 4 ) 
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where 0 is the NG boson for symmetry breaking, which is to be absorbed 

into the longitudinal mode of massless gauge boson A^. C" is gauge-equivalent to the 
Lagrangian C. Actually, if we take the unitary gauge {6 = —0), the Lagrangian C 
reduces to C. It is noted that in contrast to C which has no gauge symmetry, C" has a 
U{l)iocai symmetry: 

f ^ -|- \fNd^6 , 

< —>• tjj' = ) (5) 

[ (t> ^ cf)' = cf) + e. 

The total Lagrangian can be obtained by adding to (^) the gauge hxing term and Faddeev- 
Popov(FP) ghost term: 


Cgf+fp = -i5B{cf[A, c, c, 5,0]) , 


( 6 ) 


and is invariant under the following BRS transformation: 


5bA^,{x) 

= dtFc{x) , 

5bB{x) 

= 0, 

5bc{x) 

= 0, 

5bKx) 

= iB{x) , 

<5i?0(a;) 



= ^c(a;)0(a; 

5Bcr{x) 

= 0, 


(7) 


where c{x),c{x) are the FP ghost helds and B{x) is the Nakanishi-Lautrap held. If 
we chose f {A, c, c, B, (j)) = F{A,(j)) + |5, FP ghost is decoupled from the system. By 
integrating out B, we can rewrite gauge hxing term Cqf as 

CaF = ~F(AAf- ( 8 ) 

Here, covariant gauge and gauge corresponds to following choice for F{A,(()), respec¬ 
tively. 

F(A, 0) = d^A^ : covariant gauge , (9) 
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F{A, 0) = d>^A^ + gauge . 

(jT- 


( 10 ) 


V 


In the covariant gauge(H), NG boson 0(x) is not decoupled except the Landau gauge(.^ = 
0). On the other hand, in the case of R^ gauge, 0(x) is completely decoupled from the 
system and therefore the system can easily be treated]^. Here we choose R^ gauge. The 
total Lagrangian is obtained by adding £gf to the Lagrangian 


C!" = C + C^, 


C = '4){i^ + J)^} + -^'ip 


aipip 


2Gv 
1 


2Gc 


+ 


2G 


V 




2^ 


£ = ^(^^+ J)P; + V + 4^ 


zaxpxp , 


( 11 ) 

( 12 ) 

(13) 

(14) 


y/N 

where and DjPJ are the tree level scalar propagator and the tree level gauge boson 
propagator, respectively, which are written as 


G(o) = 


-Gc 


T>f°) = 


1 IIU ' 


jT 

' liv “I 1 " 

I '^V 

pT ^ r, - 

^ LLV — 






pL ^ 


(15) 

(16) 
(17) 


3 Vacuum Polarization 


In this section, we calculate the vacuum polarization tensor in the Euclidean space. 
By making use of the gauge invariant Pauli-Villars regularization, it is shown that the 
1-loop vacuum polarization tensor has the following form||14||: 


n^4p) = n(p)pj^, 

where n(p) can be written as 

^ -2(trl)r(2-g/2) fi (x - xV 

’ (4ir)0/2 Jo [F + (x- a:2)p2}2-°G 

^ (trl)r(2-£>/2) ,o ,D_4p (o_Dol _A\ 

3(4ir)D/2 P t 2’ ’2’ 4JV ’ 


(18) 


(19) 

( 20 ) 
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with F{a, j3,'y, z) being the hypergeometric function. Since the source J is inhnitesimal, 
we can consider J as 2J/p 1 for non-zero momentum p ^ 0. Using the mathematical 
identity for the hypergeometric function, 


—^(-^) “F(a,a-7 + l,a-/3 + l;^ 
r(/3)r(7-a) 

rTTr?- ^ F{P, P —'y + 1, P — a + 1; z ^) , ( 21 ) 

r(a)r(7-/l) 

we can expand 11 (p) in terms of 2J/p: 

n(p) = foP^-^ + f2P^-^J^ + fDp-^J^ + hP^-^J^ + 0(^ (22) 


/o 

/2 


U 

fo 


-2(trl)r(2 - D/2) 


(47r)-D/2 

2{A-D){D-1) 


B{D/2,D/2) , 


fo ) 


{D-2) 

2{Q - D){D - 3){D - 1) 
{D-2) 

2^0 

D{D - 2)B{D/2,D/2) ' 


fo , 


(23) 

(24) 

(25) 

(26) 


Doing almost the same calculation as the one in 11^,^ (p), the 1-loop vacuum polarization 
for scalar propagator no-(p) is obtained as 

n< 7 (p) = -Co + ci [ dx{.P -F (x - , (27) 

Jo 


where cq and ci are dehned by 

2 (trl)A ^-2 ^ _ 2(trl)r(2-D/2)(D-1) 

^ (4ir)»/2r(D/2)(D-2) ’ (4ir)o/2(B _ 2) 


(28) 


For p = 0, 


ncr(O) — —Co + Cl ^ . 


(29) 


For p 7 ^ 0, IIo- can be expanded under the condition 2J/p 1: 


n^(p) = -Co hop^ ^ + /i 2 P^ + hop + h^p^ , (30) 


6 



{trl)T{D/2)T{2 - D/2){D - 1) 


(47r)^/22^-2(L) - 2) 

/i2 = 2,[D — l)/io , 

/i4 = 2(D-3)(D-l)/io , 

_ {tTl)T{2-D/2)T{D/2){D-l) 

IT'D — 


/in = ^, (31) 

(32) 

(33) 

(47r)^/2r(Zl)Z12i-^(Zl - 2) ' 

When 2J/p > 1, Haip) should be expanded in terms of p/2J. However it has already 
been shown in the Thirring model that the contribution from the region 2J/p > 1 is 
sufficiently small and therefore does not affect their conclusion H. Since essentially the 
same analysis can be applied to our case, we use (^OD as long as p ^ 0. 


4 Generating Functional 


The generating functional W[J] is given by 

^w[j] ^ J x>iljViljVaVAe 

First we integrate out the fermion helds: 


iS 


,iW[J\ = f VaVABet ^ 4 


X exp 


Vn Vn 


Seff is the effective action which includes fermion loop effects: 


iSeffla, A] 

= i {p!G“V + iAi'iD-^A"'^ 



(35) 


(36) 

(37) 
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( 38 ) 


where Ga{p) and are the scalar propagator the gauge propagator, respectively, at 

the leading order of the 1/N expansion: 




(39) 

(40) 


After integrating out the fermion held, we have an inhnite number of non-local vertices 
in eg. (|38|) . W[J] is given by calculating all vacuum graphs which are made by connecting 
these non-local vertices in the 1/N expansion. 

First we restrict ourselves to the leading order of 1/N expansion, 0{N). Since the 
source J is inhnitesimal, the higher order terms with respect to J can be omitted in the 
calculation of hF[^]. We expand W[J] with respect to J and we will omit the terms higher 
than J^. Therefore we consider the order of various graphs in order to estimate which 
graphs should be included in W[J] up to 0{J^). For example, we estimate the graphs 
Fig.1-2 which are contained in the order parameter {(p{J) = j): 




Fig.l 


Fig.2 


We have 



(41) 


This graph is of order 0{J). From (^T]) , we also note that 1/N leading tadpole diagram 
Fig.2 is of order 0{J). 






Let us consider the n-point vertex shown in Fig.3 {n > 3). The case n = 2 in Fig.3 is 
the vacuum polarization and has already been absorbed into G„. This vertex is obtained 


by differentiating (^) with respect to source J: 



Fig. 3 


d 


n—1 _ 


(n — 1)! \idJ^ 


NcqJ- 


Nci 
D - ] 


J 


D-l 


-Nc^{D - 1)(D - 2) ■ ■ ■ (D - n + 1) 


— l)!(Zi) — 1) 


J" 




(42) 

(43) 


If n > D, (1^ ) has negative power in J and diverges when J = 0. This is the consequence 
of existence of infrared divergence of Fig.3 in the case of J = 0. By using this vertex, we 
find that there exists a graph contributing to W[J] at Since Fig.4 includes Fig.3 

(0(J^“"')) and n tadpoles (0(J"')), Fig.4 becomes of order 

Fig.4 ~ O(J^) . (44) 


It will be shown below that a graph which contains more than two non-local vertices 
of Fig.3 does not contribute to hF[T]. It is clear that the tadpole diagram Fig.5 is of order 
with D > 2, and hence the order of this graph is higher than the simple tadpole 
Fig.2. Let us consider a graph which contains more than two vertices in Fig.3. There is at 
least one tadpole Fig.5. If we convert this tadpole Fig.5 into a simple tadpole Fig.2, the 
order of the graph becomes lower, and finally we can reduce the graph to Fig.4. Therefore 
we conclude that the order of such a graph as contains more than two vertices in Fig.3 is 
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higher than that of Fig.4 which is 0{J^), and such a graph does not contribute to this 
order. We can easily found here all graphs Fig.6-8 contributing to W[J] at 0{J^) at the 
leading order of 1/A^ expansion; 



Fig.6 Fig.7 Fig.8 


These graphs are easily calculated by using PT|) and (^). 

Similar arguments can be applied to the next-to-leading order of 1/iV expansion. In 
the next-to-leading order, we have to consider the following non-local vertices in Fig.9-14 
similar to Fig.3. 



Fig. 12 


Fig. 13 


Fig. 14 
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Fig.9 can be calculated by differentiating the scalar vacuum polarization 


Fig.9 = -Tr 
^ 2 




d 


-Q f 

Jp 




d 




(45) 

(46) 


where 12 denotes the space-time volume, 12 — / d^x. Ga{p] and no(p) can be expanded 
in terms of J. Then Fig.9 is given by 

-If 1 h2p^~^J‘^ 1 


Fig.9 = 


1 - 


2 Jp - Co +hop^-"^ \ - Co + hop ^-2 j 

X {2h2P^-^J + Dhnp-^J^-^} + 0{J^+\ J3) , 

= -BiJ-B2J^-^ + 0{J^) , 

~ 0(J), 


where 


Si = 

B2 = 


4a{D - l)A^-2 
(D-2) 
aDho 


1 - — log(l + gs 
9s 


{D — 2)ho 


log(l + gs) , 


g — cqGs , gs — 


hoA^-2 
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a = 


(47) 

(48) 

(49) 

(50) 

(51) 

(52) 


Co yl-gj^^~ (47r)«/2F(Zl/2) 

Fig.10,11,12,13, and generally n{n > 3)-point vertex(Fig.l4), can be obtained by the same 
procedure: 

2 


Fig-10 + ^Fig.ll = 


G. 


Fig.12 + Fig.13 = 


2 ■ 3! 


Tr 


d ^ 

3 




(I/)’ 


~ 0 ( 1 ) , 

0(J^-3) 


Fig.14 = 


2-n!' 


Tr 


G„ 


^ d " 
idJ. 




0(J 


D—n\ 


(53) 

(54) 

(55) 


Thus n{> 3)-point vertex Fig.l4 as well as Fig.3 is of order and this fact leads to 

the same argument as at the leading order. Note that there appear several tadpoles(0( J)) 
Fig.15-17 in the next-to-leading order: 
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Fig. 15 


Fig.16 


Fig. 17 


By using almost the same discussion as in the leading order, we hnd the following graphs 
will contribute to generating functional in the next-to-leading order. The above tadpoles 
Fig.15-17 and the following graphs Fig.18-24 can be easily calculated by use of (|4T|), (^D]), 



(HD. 



Fig.22 


Fig.23 


Fig.24 


Now we include the gauge interaction, which can be performed by almost the same 
procedure as before except that gauge tadpoles do not exist. Here we denote only 
graphs(Fig.25-31) contributing to the generating functional: 
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Fig.28 



Fig.29 


Fig.30 


Fig.31 


In order to calculate these graphs, we need Fig.32 which is obtained by differentiating the 
vacuum polarization tensor: 



Fig.32 


Namely, 


Fig.32 = -AiJ - + O(J^) , 


(56) 


^1 


2a(B - 1)/2A^-2 

(D - 2)/o 


1 - — log(l + gv) , 
gv 


(57) 
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( 58 ) 


^ _ aD{D - l)fD , ^ 

" {D-2)fo +9v), 

where gy = —/oA^“^Gy is the dimensionless Thirring coupling constant. 

In the end of lengthy calculation, we found all graphs(Fig.6-8 and Fig.18-31) contribut¬ 
ing to the generating functional W[J] in the next-to-leading order of 1/iV expansion. By 
calculating all these graphs, we hnally obtain W[J] as 


_Vh[J] = (Fig.6) + (Fig.7) + (Fig.8) + (Fig.l8) + --- + (Fig.31) 

= -KJ"^ + —PJ^ 

2 ’ 


where K, P are given by 

K = NKo 
Co 


Ki , 


K,= 
Ki = 


i-g 

Ai + Bi 


P = iVFo + A , 


Pn = 


-Cl 


D 


D - 1 
Pi = 1—^2 


B 2 + 


ciD{Ai -|- Pi) / g 


D 


Co(P-l) \^-g, 


g, 


(59) 

(60) 

(61) 

(62) 

(63) 

(64) 

(65) 

( 66 ) 


5 Effective Potential by Inversion 

The order parameter for the dynamical symmetry breaking, (p(J), is given by differ¬ 
entiating (0) with respect to source J: 

1 dW r. , 

■e(J) = = jjgj = Kj + . (67) 

We here made use of inversion methodto obtain the effective potential instead of 


taking a Legendre transformation. By inverting (|^) , we obtain 


FT-V - K-^PJ^~^ 

(68) 

K-^cp - PK-^ip^-^ 

(69) 

K~^ip + . 

(70) 
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K ^ is expanded in 1/A^ as 


From the relation 


K-^ = 


NKa 


J = 


Ki 

NKn 


dV 
^ ’ 


(71) 


(72) 


the effective potential of Zl(2 < 77 < 4)-dimensional Thirring-Gross-Neven model can be 
obtained in the next-to-leading order of 1/iV expansion: 


r = i/i'-V" + ■ 


(73) 


Dynamical symmetry breaking occnrs if the eqnation J = ^ ^ has a non-zero solntion, 

(p ^ 0. It is realized when K~^ < 0, and the critical snrface is given by K~^ = 0. Note that 
we calcnlate all graphs throngh vacnnm polarization. Therefore onr effective potential and 
all onr conclnsion has no gange dependence [|I^. 


6 Critical Line 

77“^ is given in ([7l|) reads 

77-' = 


NKo 
^-9 
Nco 
1 ( 


iVcn 


(l - ) 

V NKoJ 

(74) 

1 ^ iVco (1-g)^ f 

(75) 

^ Nco \ ' 

(76) 

given by the equation 77“' 

= 0 which dehnes a 


snrface((7c, fi'Vc) "^c) in the parameter space {g,gv,N): 


0 — 1 ~ 5'c + 


P-1) 

2Nr 


1 , / 

1 - — log(l 

9 sc 


{D-lf{A-D) 
2N,{D - 2) 


1-log(l + 9vc) 

9vc 


77) 

where gsc was dehned in (^) . We can see from (O) that (yfc ~ 1 for large N^, and therefore 
we can approximate this eqnation as 


^ 1 + 


77-1 (77-1)2(4-77) 

2Nc 2Nc{D - 2 ) 


1-log(l gvc) 

9vc 


+ 0 ^logiVj . (78) 
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Eq.([75D works well when N is sufficiently large. We will investigate the structure of the 
critical surface using eq.( 


Before investigating (^) in detail, we consider the pure Gross-Neveu model. Taking 


gvc = 0 in the eq.([78D, the critical line of the dynamical symmetry breaking in the Gross- 
Neveu model is given by 


It is remarkable that this critical line is exactly the same line obtained by using SD 
equation up to the next-to-leading order of 1/iV expansion ||I6|. 

Let us return to ([78| ) to look into the critical surface of Thirring-Gross-Neveu model. 
We consider here D = 3 for convenience. Then the critical line is written as 


9c 


. 1 ~ ‘^H{gvc) 

Nc 


( 80 ) 


where H{z) is dehned by 

H{z) = 1--\og{l + z) . (81) 

H{z) is monotonically increasing function and H{0) = 0,H{oo) = 1. Since 0 < iL < 1 
and > 1, we hud Qc < 2 in (|80[) . The critical line in the {gv,N) plane is given in 
FIG.A-1 and FIG.A-2 for various values of g. 

We here divide the region 0 < 5 ^ < 2 into three regions, (l)^' > 2, (2)0 < g < 1(FIG.A- 
1) and (3)1 < 5 ^ < 2(FIG.A-2), where g = 1 is the value of critical coupling constant at 
the leading order of the Gross-Neveu model. The forms of critical line in these region are 
different from each other as follows: 


( 1 ) <7 > 2 

In this case, the chiral symmetry is always broken irrespectively of gv and N. 
{2)l<g<2 (FIG.A-1) 

In this region, the critical line (|5(]|) reads 

5'c - 1 fi-C - 1 


(82) 


For g > 1, gvc can take a value which satisfy the equation 2H{gvc) < 1- The critical line 
exists only in the region N < and Nc decreases as gvc increases as shown in FIG.A-1. 
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The symmetry is dynamically broken when N > 

(3) 0 < c/ < 1 (FIG.A-2) 

In this case, Nc increases according to increasing gvc and there exists a critical value 
of TV,: 


Nc = 


(gvc) — 1 
1 - 5'c 


< N^{gc,gvc = oo) = 


(83) 


The dynamical symmetry breaking does not occur when N > Nc{gc,oo). This critical 


line is quite similar to the critical line obtained in the Thirring model |T3 . Thirring 
interaction strongly affect the system in this region, 0 < 5 ^ < 1. Thus, the property of the 
critical line is drastically changed at g = 1. 

This can be understood naturally as follows. 5 ^ = 1 is the critical coupling constant of 
the Gross-Neveu model at the leading order of 1/A^ expansion, and Thirring interaction 
does not strongly affect the leading order. Therefore, for g > 1, Gross-Neveu interaction is 
dominant. The system behaves like the Gross-Neveu model, and the symmetry is broken 
almost all the region. As N increases, the effect of 1/A^ next-to-leading order becomes 
small, and the system is dominated by the tadpole of the leading order. Therefore the 
dynamical symmetry breaking occur in large N region for g > 1. 

On the other hand, the critical behavior is quite different in 0 < ^ < 1 . In this region, 
the Thirring interaction strongly affects the system, and the system (or the critical line) 
becomes Thirring-like. There exists a critical value for N, Nc{gc,gvc = cxo) < cx). The 
dynamical symmetry breaking does not occur in the large N region, N > Nc{gvc = 00 , gc), 
even if gy —> 00 . The symmetry is broken in the small N region. This is in contrast to 
the case 1 < 5 ^ < 2 where the system behaves like Gross-Neveu model. 

Thus the critical behavior drastically changes at g = 1. The system behaves like 
Thirring model in < 1, whereas it does like Gross-Neveu model in > 1. 


7 Conclusion and Discussion 

In this paper, we have constructed the effective potential for the order parameter of the 
chiral symmetry, the fermion condensate {'ijj'ip) in the hybrid model of the Gross Neveu 
model and the Thirring model in D{2 < D < 4) dimensions. From this potential, we 
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have shown existence of the chiral symmetry breaking and explicitly obtained the critical 
surface. 

We have analysed our model using the inversion method in 1/A^ expansion. In this 
model, existence of Gross-Neveu interaction make our analysis complicated. So many 
diagrams contribute to the generating functional. But we found we do not have to cal¬ 
culate all the diagrams. We have expanded the generating functional W[J] in terms of 
the inhnitesimal source J, and estimate it on the order 0{J^). In section^, we investi¬ 
gated systematically the order of diagrams with respect to source J, and determine which 
diagram contributes to the generating functional W[J] or effective potential V{{'ip'ip)). Re¬ 
sultant effective potential has no gauge dependence, because we calculated all diagrams 
through vacuum polarization. 

From this effective potential, we obtained explicitly the critical surface which separates 
symmetric phase and broken phase of the chiral symmetry. Especially, if we consider the 
special case gv = ^ (Gross-Neveu model), our result ([79|) reproduces exactly the same 
result obtained by using SD equation||l6|. 

Kim et al.[^ have already studied Thirring Gross-Neveu model by using SD equation 
and they found critical surface. Although the original Lagrangian(2) has no gauge sym¬ 
metry, they pretended the auxiliary massive vector held as really a gauge held without 
introducing hidden local symmetry. Furthermore, the SD equation becomes coupled in¬ 
tegral equations for A{p) and B{p) in the full fermion propagator S{p)~^ =']^A{p) —B{p), 
which are difficult to solve. Therefore additional approximations (or assumptions) be¬ 
sides 1/N expansion were made in their analysis. They assumed that the form of A{p) 
can be determine by using perturbation, and mass function B{p) can be regarded as a 
constant. As a result of such an approximation, the ehects of the next-to-leading order 
of the Gross-Neveu interaction were not included. 

On the other hand, we succeeded in obtaining the ehective potential and the critical 
line without using additional approximation besides 1/N expansion. 
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